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ABSTRACT 

We examine galaxy formation in a cosmological AMR simulation, which includes two 
high resolution boxes, one centered on a 3 x 10 14 M cluster, and one centered on a 
void. We examine the evolution of 611 massive (M* > 1O 1O M0) galaxies. We find that 
the fraction of the final stellar mass which is accreted from other galaxies is between 
15 and 40% and increases with stellar mass. The accreted fraction does not depend 
strongly on environment at a given stellar mass, but the galaxies in groups and cluster 
environments are older and underwent mergers earlier than galaxies in lower density 
environments. On average, the accreted stars are ~ 2.5 Gyrs older, and ~ 0.15 dex 
more metal poor than the stars formed in-situ. Accreted stellar material typically lies 
on the outskirts of galaxies; the average half-light radius of the accreted stars is 2.6 
times larger than that of the in-situ stars. This leads to radial gradients in age and 
metallicity for massive galaxies, in qualitative agreement with observations. Massive 
galaxies grow by mergers at a rate of approximately 2.6%/Gyr -1 . These mergers have 
a median (mass- weighted) mass ratio less than 0.26 ± 0.21, with an absolute lower 
limit of 0.20, for galaxies with M* ~ 10 12 M . This suggests that major mergers do 
not dominate in the accretion history of massive galaxies. Al l of these re sults agree 
qualitatively with results from SPH simulations by lOser et al.l (|2010l |2012[ ). 
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1 INTRODUCTION 

In a ACDM cosmology, galaxy growth is hierarchical; large 
galaxies form by accreting smaller systems of gas, stars, and 
dark matter. This model has been extensively tested by sim- 
ulations and observations. However, it is not known whether 
stars or gas dominate the accretion history of galaxies, and 
what effects this accretion has on the observable proper- 
ties of the galaxies. With modern hydrodynamical simula- 
tions, we can directly measure what fraction of the stars in a 
massive galaxy have been formed in-situ and what fraction 
have been made in other galaxies which were subsequently 
accreted by the parent galaxy. Similarly, we can map out 
when accretion occurs, what kind of accretion events domi- 
nate, and how the accretion affects the observable properties 
of galaxies. In this work, we explore these effects for a set of 
massive galaxies from a hydrodynamical cosmological simu- 
lation. 

lOser et all (|2010h (hereafter OslO) present a useful 
framework for studying the hierarchical build-up of galax- 
ies. Based on hydrodynamical simulations of early- type 
galaxy formation, they suggest that massive galaxy forma- 
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tion can be divided into two phases: an early, rapid in-situ 
star formation period followed by a late merger-dominated 
period. The former period is similar t o the picture pre- 
viously called "monolithic collapse" (jEggen et al. I Il962l : 
iLarsonl Il975l : ICarlberd 1 19841 ) and currently labelled "cold 
flow- driven" star forma tion (e.g. iDekel &; Bir nboim 2006; 
I Dekel et al.1 120091 : lElmegreen et al.1 120091 : ICeverino et~ai1 
while the latter phas e is similar to prior pic- 
tures of galactic cannibalism ( Ostriker & Hausman 1977; 
White fe Reeslll978l : ISearle fe Ziimlll978l : iKauffmann et al.1 
19931 ). In this scheme, individual stars can be classified ac- 



cording to whether they formed within the virial radius 
of the final system (in-situ) or outside of it (accreted). 
OslO show that this distinction between stars is unambigu- 
ous; there is a large spatial separation between the birth 
places of stars formed in-situ and those added in accre- 
tion events. They further show that massive galaxies are 
dominated by accreted stars, which are typically older than 
stars formed in-situ. In this work, we find that the ac- 
creted stars in a galaxy are older, more metal-poor, and 
found at larger radii in the final galaxy than stars formed 
in-situ. The metallicity trends are co nsistent with observa- 
tions of radial trends in colour (e.g. de yaucouleursl Il96ll : 



iTortora et all l20ld : iTal fe van Dokkuml 1201 if ), metallicity 
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(Spinrad fc Tavlorl Il97ll: 


Faber et alJ Il977l; iDavies et al. 


1993: Brourfi et all 120071: iRawle et alJ 12008: Soolaor et al. 


2010l:lKuntschner et al.ll20: 


LOl:ICoccato et al.ll201lh and glob- 


ular cluster metallicities ( 


Forbes et alJ l201ll: lArnold et al. 



2011) in early- type galaxies. 

The differences between accreted and in-situ stars are 
also seen in disc galaxies. For example, in the Milky Way 
stellar halo, there has long been observational evidence for at 
least two stellar populations. The inner h alo consists of stars 
with high a— element abundances (e.g. iNissen &; Schuster 



2010) and rotates with the galactic disk (De ason et al.l 
201l|), while the outer halo has low a-element abundances 
and is not rotating, suggesting it was created by accretion 
of satellite galaxies. Simulations of disc galaxies show that 
the differences in dynamics and chemical abundances in the 
stellar halo are dues to the differences in the accreted and 
in-situ stellar populations (iBrook et al.l 2004 MAbadi et al.l 
120061 : IZolotov et ai1l2009l . l20ld : lFont et al.ll201lh . " 

The two phase model of galaxy formation (i.e., early, 
in-situ star formation and late accretion) is also useful 
in explaining the evolution of comp a ct, massive ellipti- 
cals at z ~ 2 dTrujillo et al.1 [20071; van Dokkum et al. 
20081: iDamianov et al.ll2009l : Ivan Dokkum fe BrammerfeOlOl : 
Oser et al.1 l2012h . Both direct profile measurements and 



velocity dispersions have shown that these systems are 
100 times denser (within 1 effec tive radius) than 



present-day ellipticals of t he same mass (Daddi et al. 



2005 



2011 



van Dokkum et al.l I2OO8J ; Ivan der Wei et al.1 I2008T 
van de Sande et al.l l201lh . However, if compared on same 
physical scale, the densities in the central portions of z ~ 
2 and present-day ellipticals are similar, suggesting that 
early-type galaxies have increased their size through mi- 



sive galaxies (iNaab et al.1 120071. 12009 


; iBezanson et al. 2009; 


Hopkins et al.l 120091; ICarrasco et al.l 


2010:bser et al.lboi2; 


Tal et al.ll2012l, but see Newman et al. 2012^. Bv using close 



companions to es timate the merger rate of compact galaxies 
at 0.4 < z < 2.5, iNewman et al.l (|2012n observe that minor 
merging may not be sufficient to explain the size evolution 
of compact galaxies from z ~ 2. However, the size evolution 
due to minor m ergers has been found in hydrod ynamical 
simulations (e.g. INaab et al.l [20091 : lOser et al.ll2012h and is 
consistent with the size growth found by observations. 

In this work, we make use of galaxy catalogs from 
cosmological si mulations done with adaptive mesh refine- 
ment (AMR) ([Cenl l2Q12h . The use of an Eulerian grid- 
based code instead of a Lagrangian particle-based code is 
a notable distincti o n b etween this work and that of OslO. 
Scannap ieco et al. I (|201lh find that different numerical hy- 
drodynamics can yield differences of a factor of two in simu - 
lated galaxy properties. However, IScannapieco et all (|201lh 
find that changes in the feedback implementation can yield 
even larger differences in the simulation results. We will ad- 
dress the differences in sub-grid physics between our work 
and that of OslO below. The AMR simulation studied here 
contains two high resolution boxes, one centered on a galaxy 
cluster and another centered on a void. These large volumes 
simulated at high resolutions allow us to study the merg- 
ing histories of more than 600 galaxies in a variety of en- 
vironments, from void to group to cluster. Following OslO, 
we have divided the stars in each galaxy into accreted and 
in-situ, which allows us to examine the star formation his- 



tories for the galaxies as a function of stellar mass and en- 
vironment. The qualitative similarities between our results 
and those of OslO, despite profound differences in numer- 
ical techniques and sub-grid physics, help substantiate the 
two-phase model of galaxy evolution for early type galaxies. 

The paper is divided into the following sections: 321 de- 
tails the simulation, the building of the merging histories 
and the tagging of stars as accreted or in-situ; 33 and 33.11 
describe the properties of the in-situ and accreted stars and 
compare our simulation results to those of OslO. 2] focuses 
on the merger histories and the types (major vs. minor) 
of mergers the galaxies in our simulation underwent and 
their effects on the final galaxy properties. For the simula- 
tion and throughout this work, we us e the following ACDM 
cosmology, consistent with WMAP-7 (Komatsu et al ] l201lh : 
Q M 0.28, Q A 0.72, Q b 0.046, a 8 0.82, n 92, 
Ho = 100/ikm s" 1 Mpc" 1 , and h = 0.7. 



2 DESCRIPTION OF SIMULATION 



This work uses the same cosmological simulations as Cen 
(|20ll l2011bh . and further descriptions of the methodol- 
ogy can be found there. The simulations are perfo rmed 
with the AMR Eulerian hydrodynamics cod e, Enzo (|Brvanl 
1 19991 : lO'Shea et ai1l2004l : Ijoung et al1l2009h . Enzo uses the 
particle-mesh technique for the dark matter N-body dynam- 
ics. Poisson's equation is solved using a fast-Fourier trans- 
form on the base grid and mult igr id relaxatio n on the the 
higher resolution subgrids ([O'Shea et al.ll2004h . Enzo solves 
the hydrodynamic s using a modified version of the piecewise 
parabolic method (jWoodward fc Colellal[l984l : ICT Shea et al. 
2004), which is formally second-order accurate and conserves 
energy, momentum and mass fluxes. For this work, two re- 
gions, a cluster and a void, were selected from a low reso- 
lution simulation of a periodic box 120 h~ 1 Mpc on a side. 
These regions were re- simulated separately at high resolu- 
tion. The cluster box is 21 x 24 x 20 /i~ 3 Mpc 3 , and the void 
box is 31 x 31 x 35 /i _3 Mpc 3 . These regions represent +1.8cr 
and — l.Ocr fluctuations, respectively. The refined regions are 
surrounded by two buffer layers, ~ l/i -1 Mpc thick, which 
connect to the underlying low-resolution periodic simulation 
box of 120/i _1 Mpc. The resolution in the refined regions is 
460 h~ 1 pc physical or better. The two refined regions are 
similar to t wo of the five regions taken from the Millennium 
Simulation (ISoringel et alj|2005h and then simulated at high 
resolution in Crain et al.1 (|2009h . The volume of the void box 
used here is 5 tim es larger than the regions simulated in 
ICrain et al.l ([2009), whil e the cluster box is comparable in 
size to the +2cr region in lCrain et al.1 (|2009h . 

In the refined regions, the dark matter particle mass 
is 1.07 x 10 8 M©, and the stellar particle mass is typi- 
cally 10 6 M©. Each stellar particle is tagged with its ini- 
tial mass, metallicity, dynamical time, and formation time. 
Additionally, the luminosity of each star particle in the 
five S loan Digital Sky Survey (SDSS) filters (jFukugita et al.1 
1996) is computed using the Galaxy Isochrone Synthesis 
Spectral Evolution Librar y (GISSEL) stellar synthesis code 
(jBruzual fe Charlotll2003h . assuming a Salpeter initial mass 
function (IMF). 

The simulat i ons include a UV background 
(Haa rdt &; Madaul Il996h , a prescription for shielding 
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from UV by neutral hydrogen (|Cen et aL 200 j), and met al- 
licity dependent radiative cooling i|Cen et all Il995h. The 
UV background radiation field is subject to self-shielding; 
we compute the local ionization and cooling/heating 
balance for each cell in the following way. A local optical 
depth approximation is adopted to crudely mimic the local 
shielding effects: each cubic cell is flagged with six hydrogen 
"optical depths" on the six faces, each equal to the product 
of the neutral hydrogen density, hydrogen ionization cross 
section and scale height. The appropriate mean from the 
six values is then calculated. Equivalent means for neutral 
helium and singly- ionized helium are also computed. Metal 
cooling is computed with ionization equilibrium including 
the effects of the redshift-dependent UV- X-ray background, 
using a code based on the standard Raymond-Smith code 
([Raymond et al.1 Il976h . The cosmological reionization in 
this model occurs at redshift z — 8 — 9. 

Star par ticles are created in cell s that satisfy the criteria 
set forth in ICen fe Ostrikerl (|l992h . Specifically, stars are 
created when the baryon overdensity, 5pb/p > 5.5 and the 
following are true: 



V • v < , 

^-cool ^ ^-dyn — 
TUb > mjeans 



3tt 



32Gp t . 



, and 



G^pl /2 (l + p b /p d )^ 



(la) 
(lb) 

(lc) 



where c s is the sound speed and pb and pd are the baryon 
and dark matter densities, respectively. 

Supernovae feedback is modeled as in lCen et al.l (|2005l ). 
Feedback thermal energy and metal-enriched eject a is dis- 
tributed into the 27 cells surrounding the star particle. 
The amount added to each cell is inversely proportional 
to the density of each cell, mimicking the process of su- 
pernovae blast waves, which channel energy into the least 
dense surrounding regions. The rate of feedback has the form 
1/tdyn [(£ - U)/tdyn] exp [—(t - U)/tdyn], where U is the for- 
mation time of a given star particle, and td yn the dynamical 
time of the cell. Accounting for both Type II and prompt 
Type I supernovae, the total amount explosion kinetic en- 
ergy per stellar mass (M*) formed (assuming a Chabrier 
IMF) is e SN M*c 2 , where e S N = 1 x 10~ 5 . The feedback pro- 
cess is entirely hydrodynamically coupled to the surround- 
ings and subject to all physical processes, such as cooling. 
The simulations do not include input from active galactic 
nuclei; its absence may contribute to the excess of massive 
galaxies. 

Several tests of the feedback process for these simula- 
tion are presented in earlier works. [Cenl (|2012h show that the 
kinematic properties of damped Lyman-alpha systems in the 
simulations agree well with observati ons, providin g an sig- 
nificant test to the feedback model. In lCenl (|2011bh we show 
that the phenomenon of cosmic downsizing can be naturally 
explained in our model with the necessary large- volume cos- 
mological setting, due to correlations between galaxy envi- 
ronment an d as t rophy sical processes that operate on gas 
physics. In ICeJ (|2011ch we show that galaxy luminosity 
functions for both UV and FIR selected galaxies can be self- 
consistent ly produced by the simulation, when reprocessing 
of optic al- UV radi ation into IR radiation is followed in de- 
tail. In[Cen| (|2011aT ). we show that our model reproduces 



a variety of observed properties of O VI A A 1032, 1038 ab- 
sorption lines at z ~ 0. In combination, these tests strongly 
indicate a range of applicability of our simulations to com- 
plex systems, from galaxies at kpc scales to Mpc scales in 
the intergalactic medium, and empirically validate our star 
formation feedback prescription, upon which all these phe- 
nomena critically depend. This validation of the simulation 
results is critical and allows us, with significant confidence, 
to perform the particular analysis presented here. 

2.1 Merger trees 

In order to identify galaxies in a given snapshot, we use 
the HOP algorit hm on stellar particles ([Eisenstein Sz HvJ 
ll999l ; ICenll2011br ). For the HOP program we choose an outer 
stellar overdensity of 6 outer = 10000, and density ratios of 
<5 'saddle = 2.5(5 Q nter and S vea k — pouter ', these ratios are rec- 
ommended by lEisenstein &; Hul (| 19991 ). We have examined 
600 groupings with M* > 5 x 1O 9 M at z — in the cluster 
simulation and have found 9 HOP groupings which contain 
more than one (usually two) visually distinguishable galax- 
ies. The central galaxy of the cluster is included in one such 
grouping. We find no such groupings in the void box. This 
represents a contamination of < 2% of our galaxy sample. 
For the remainder of the paper, we will exclude the 9 group- 
ings which have not been adequately separated by the the 
HOP algorithm at z — 0. 

In order to study the accretion histories of the galaxies, 
we build merger trees for the simulations. For the cluster 
box, we use 38 snapshots: one every Az = 0.05 between 
z = and 1.35. Beyond z — 1.35, we have snapshots at 
z = 1.5, 1.6, 1.75, 1.9, 2.0, 2.2, 2.5, 2.8, 3.1, and 4. For the 
void box, we only use 14 redshift slices at z = 0, 0.05, 0.15, 
0.2, 0.4, 0.5, 0.6, 0.8, 1.0, 1.6, 1.9, 2.5, 3.1 and 4. The merger 
trees are built by looking at the overlap of the stellar parti- 
cles in galaxies in sequential redshift slices. For example, if 
most of the stars in galaxy B at z — 0.05 are later in galaxy 
A at z = 0, we label galaxy B a parent of galaxy A. Galaxy 
A may have many parents at z — 0.05; we call the main pro- 
genitor (the same galaxy at an earlier time) the galaxy that 
contributes the most mass to the child galaxy. The remain- 
der of the parent galaxies are galaxies which have merged 
into the main galaxy during the current redshift step. This 
process can be repeated for all the galaxies in all the red- 
shift slices. Each galaxy has a progenitor and possibly a set 
of galaxies which have just merged into it. Stepping back 
through the main progenitors gives a full history for a given 
galaxy. 

We verify this method of grouping by using every other 
redshift slice and checking that the progenitors do not 
change. Ninety-two percent (91%) of the galaxy merger trees 
in the cluster (void) box are the same to z — 1 if we only use 
every other redshift slice. The galaxies with modified histo- 
ries are typically below our mass threshold. We also check 
that the 200 most-central particles in a galaxy have signifi- 
cant overlap with the 200 most-central particles in its main 
progenitor. This is true for more than 98% of the galaxies 
with masses above 5 x 10 9 for the cluster box and 95% of 
the galaxies in the void box. The majority of galaxies which 
fail this test cannot be traced back to z > 2.0 and will not 
be used in our analysis. 

After identifying merger histories for each galaxy, we 
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can label the stars in every z — galaxy as accreted or 
in-situ. As above, we proceed in a stepwise manner toward 
high redshift snapshots. For each stellar particle in a galaxy 
at a given redshift, we examine where it was in the previous 
redshift slices. Star particles which were in the main progen- 
itor galaxy are labelled as in-situ. Stellar particles in any of 
the other parent galaxies are labelled as accreted, and the 
accretion redshift is noted. The remainder of stellar parti- 
cles were either formed during the current redshift step, and 
labeled as in-situ, or they formed earlier and are labeled as 
accreted. These stellar particles are due to accretion from 
unresolved galaxies, as well as uncertainty in the grouping 
algorithm. If a star particle was born in the main progen- 
itor galaxy, but later re- accreted, we consider it an in-situ 
star particle. Such particles contribute between 2 and 5 per 
cent to the total galaxy stellar mass. We then examine the 
stars of the main progenitor (those which were initially la- 
beled in-situ) in the same way and repeat the process for all 
higher redshift slices. This allows us to distinguish between 
accreted and in-situ stellar material to our highest redshift 
snapshot at z = 4. 

One difficulty in this method is ensuring that none of 
the stars we consider in-situ were actually formed outside 
the galaxy and accreted between snapshots. OslO find a 
large separation in space (3 orders of magnitude in dis- 
tance from the centre of the main progenitor galaxy) be- 
tween in-situ and accreted particles. However, their Figure 
7 also shows the birth places of a population of accreted 
stars which formed between 0.1 r v iriai and 1.0 rviriai- If our 
redshift snapshots are not sufficiently close together, we will 
count these stellar particles as in-situ, instead of accreted. 
We have compared the dynamical time of a virialized halo 
(tdyn ~ l/VG 200p C rit ~ 1/ (10H (z))) , to the time between 
redshift slices. For the cluster, the dynamical time of the halo 
is larger than the time between redshift slices by at least a 
factor of 2 for z ^ 1.35. Beyond z — 1.35, the halo dynami- 
cal time is comparable to the density of redshift snapshots. 
Therefore, our calculations of the accreted fraction may not 
include stars formed during infall to the final galaxy, and 
should be considered lower limits. 



2.2 Galaxy population 

Before examining the merger histories of the simulated 
galaxies, we have exam i ned th e overall properties of the 
galaxies at z = 0. ICe3 <|2011bh have shown that the lu- 
minosity function for galaxie s in the cluster and void agree 
well with SDSS observations ([Blanton et al.ll2003l ), except at 
the highest luminosities, where the simulation produces too 
many luminous galaxies. This is due to insufficient feedback; 
the simulations do not include AGN feedback or stellar feed- 
back from Type I supernovae. Add itionally, the highest mass 
galax ies suffer from overmerging ( White] 1 19761 : | Moore et aU 
1996). These highest luminosity bins can be corrected by 
taking into account AGN feedback and decreasing the stel- 
lar mass by a factor of 1/( 1 + (M ha io/l x 10 13 M ) 2/3 ) (jCenl 
l2011bl : ICroton et al J 12000 ). 

Figure [1] shows the uncorrected mass function of galax- 
ies in each region at z = 0, compared to the mass func- 
tion of SDSS galaxies (|Li fe White! l2009h . We have only 
included galaxies with at least 100 dark matter particles 
(Mhaio > 1O 1O M ). The average mass function for the simu- 



lation is a weighted su m of the clus ter and void box. Using 
the same simulations, ICenl (|2011bh suggest a weighting of 
6 : 1 for the void and cluster box, based on matching the 
simulated star formati on rate hi s tory t o observations. The 
weighting scheme of ICrain et al.l (|2009h , yields a weighting 
of 2 : 3 for the void and cluster box, since the cluster box is 
63% smaller than the void box, but it's also a rarer density 
perturbation. For the remainder of this work, we weight the 
galaxies in the void and cluster box equally, unless otherwise 
noted. 

There is a clear deficit of galaxies with masses below 
~ 3 x 10 10 M©. In fact, if we include galaxies with fewer 
dark matter particles, both mass functions turn up again 
and have second low mass peaks at M* ~ 10 8 M©, how- 
ever, these groupings are most likely due to resolution ef- 
fects. The simulated mass functions bracket the observations 
well in the range 10.5 < logM*/M© < 11, but as explained 
above, the simulation overproduces massive galaxies. The 
overproduction of massive galaxies is a typical problem for 
hydrodynamical simulations and is due to insufficient feed- 
back to suppress star formation and to overmerging among 
close pairs of galaxies. This is a serious limitation of our 
work, and the numerical values we derive from the simula- 
tions will certainly change with higher resolution, improved 
feedback models, and lower star formation rates. 

Since improved feedback will reduce the stellar masses 
of all massive galaxies, it is not obvious if it will decrease 
number of accreted or in-situ stars more. The high levels of 
star formation at z = (see £13. 1[) suggests that most of the 
excess stars are formed in-situ, but only improved simula- 
tions will give a definitive answer. Nonetheless, we are confi- 
dent in our results for two main reasons. The first is that the 
differences between accreted and in-situ stars do not depend 
strongly on mass (see Figures [5HZJ) ? and are present even in 
the lowest mass bin where the mass function of galaxies is 
reasonable. Secondly, although improved feedback will cer- 
tainly modify the mass function, it will not substantially al- 
ter the rank-ordering of galaxies; the most massive galaxies 
in our simulation will still be the most massive. Therefore, 
we expect that trends with stellar mass will have different 
slopes, but the sign of those slopes should be insensitive to 
the feedback model. In the next sections, we will only exam- 
ine the properties of galaxies with at least 100 dark matter 
particles and stellar masses larger than 3 x 1O 1O M , but we 
continue to include the most massive galaxies despite their 
overabundance. 



2.2.1 Galaxy formation efficiency 

Simulations from cosm ological initial conditions are kn own 
to overproduce stars ([Oser et al.ll2010l : iGuo et al. I l2010l and 
references therein). This is true for the simulations in this 
work as well, and can be easily shown by comparing the mass 
ratio of dark matter and stars to observations. The ratios of 
stellar mass to dark matter mass at z = are 1.2 x 10 ~~ 2 
and 6.0 x 10 ~ 2 for the void and cluster boxes, respectively. 
The two weighting schemes discussed above give an average 
ratio in the range 1.9to4.1xl0 -2 for the entire simulation. 
The obser ved ratio in the local universe is (1.2 ± 0.3) x 
10" 2 (e.g. iFukudta fe Peebles! 120041 : iGallazzi et al] l20Q8h . 
The void box has a star formation efficiency which agrees 
well with the observed average. However, observations show 
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Figure 1. The z = mass function of the cluster and void box 
shown in red and blue, respectively. The mass function from SDSS 
is shown by the black dashed line dLi W hite 200^). The mass 
functions of the simulated boxed bracket the observations for 
10 10 M© < M* < 10 1:l Mq. The simulations overproduce galaxies 
more massive than 10 11 Mq. 



that underden sities in the universe are also under- luminous 
(Peebles 2001). Therefore, if the star formation and feedback 
in the simulation were accurate, the void box would have a 
lower star to dark matter ratio than the average. Thus, stars 
are overproduced in both the cluster and the void box, with 
the overall star formation rate too high by a factor of two 
to four. 

We can examine the overproduction of stars in the 
cluster box in more detail by computing the ratio /* = 
M*/Mdm x (QdmM) for both simulated and observed 
galaxies and their host dark matter halos. In the simulations 
used here, the central cluster has a virial radius (using the 
radius inside which (p) — 200p C rit) of 1.3 /i _1 Mpc. The total 
dark matter mass interior to this radius is 3.0 x 10 14 Mq, and 
the total stellar mass interior to this radius is 3.0 x 10 13 M©. 
Therefore, the efficiency of star formation in the cluster 
is /* = 0.60, or 60% of the baryons in the cluster have 
been turned into stars by 2; = 0. For a cluster of this 
mass, the expected star formation efficiency from weak lens- 
ing and halo occupati on distribution (H OD) methods is 
10% < /* < 15% (|Leauthaud et alJ l2Q12h . This range is 
slightly above the values found by matching simu lated DM 
halos to the observed galaxy mass functions (Guo et al.l 
l20ld : iBehroozi etHI l20ld : iMoster et aD l20ld ). Thus, the 
star formation in the cluster is over-efficient by roughly a 
factor of four. The simulation star formation efficiency is in 
part too high due to the lack of sufficient feedback from su- 
pernovae and AGN. If we correct the stellar masses to take 
this into account, as in the previous section, the efficiency 
decreases to ~ 30%. 



2.2.2 Mean stellar properties 

The top panel in Figure [2] shows the relation between 
galaxy stellar mass and SDSS r— band luminosity- weighted 
age. The simulated trend is compared the observed mass- 
age relation from SDSS (thick gray line) (Gall azzi et al.l 
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Figure 2. The mean (luminosity- weighted) stellar ages for all 
the stars (top), and the in-situ (middle) and accreted (bottom) 
stars as a function of total galaxy stellar mass. The large black 
triangles are the medians and inter-quart ile ranges in 3 mass bins. 
The dashed line is fit the medians. The gray solid line in the 
top panel is the me dian observed mass-age relation from SDSS 
(Gall azzi etafl 12005). The galaxies in the simulation are young 
compared to observed galaxies. Nonetheless, the accreted stars 
are uniformly older than the in-situ stars. 



l2005h . Although there is little trend with mass, the youngest 
galaxies are among the least massive. The lack of trend 
at high masses is co nsistent with observations from SDSS 
(|Gallazzi et al.1 2005) , in which the steepest portion of the 
stellar mass-stellar age trend occurs in the mass range 
3 x 10 9 < M* < 3 x 10 10 . However, the median ages for 
the galaxies in our s ample are several G yrs younger than 
observed galaxy ages ([Gallazzi et alJl2005h . This is partially 
due to an excess of star formation at late times (see Figure 
[5} which is related to the insufficient resolution and feedback 
modeling in the simulation. Figure[2]also shows that the me- 
dian age of galaxies is dependent on environment, with the 
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Figure 3. The mean (luminosity- weighted) stellar metallicities 
for all the stars (top), and the stars divided into in-situ (middle) 
and accreted (bottom) stars as a function of total galaxy stellar 
mass. The large black triangles are the medians and inter-quart ile 
ranges in 3 mass bins. The dashed line is fit the medians. The gray 
solid line in the top panel is the median o bserved mass-met allicity 
relation from SDSS (Gallazzi et al. 2005). The galaxies in the sim- 
ulation show no mass metallicity relation, although the maximum 
met allicity is consistent with observations. The symbols are as in 
Figured 



oldest galaxies found in the cluster box, while none of the 
void box galaxies are older than ~ 7 Gyr. Thus the trends 
of galaxy age with stellar mass and local density are broadly 
consistent with the observational picture of downsizing and 
the quenc hing of star formation in cluster environments (see 
ICe3l2011bL for a study of downsizing in this simulation). 

It is expected that the trends in metallicity will be sim- 
ilar to those for stellar ages. As with the stellar ages, we find 
that stellar mass and stellar metallicity are not strongly cor- 
related in the simulations (see Figure [3] top panel) . There is 
a stronger dependence on metallicity with environment; the 
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Figure 4. The fraction of accreted material (M acc /M*) as a func- 
tion of final (z = 0) stellar mass. In the top panel, the blue(red) 
points are for the void (cluster) box. The large black triangles are 
the medians and inter-quart ile ranges in 3 mass bins. The line 
is fit through the medians. The lower panel shows the same for 
the galaxies simulated in OslO (data from L. Oser). The accreted 
fraction found in this work is significantly smaller than that in 
OslO for very massive galaxies. 



void box galaxies have metallicities around Z « 1 Zq , while 
the cluster galaxy metallicities are centered at Z « 1.5 Zq 
with very little scatter. This value of metallicity is con- 
sistent with the median metallicity f or the most massive 
(M > 10 11 M ) nearby SDSS galaxies (jGallazzi et alJl2005h 
(grey solid line). 

The excess of young stars and overabundance of metals 
suggests again that the si mulations suffer f rom o ver-cooling 
and insufficient feedback. IWiersma et all (|2Q1 lh find that 
changes in the feedback prescription can alter the stellar 
metallicity by factors of 4. Their simulations suggest that 
the inclusion of AGN feedback and momentum-driven super- 
novae winds can bring the stellar metallicity into agreement 
with observations. 



3 IN-SITU VS. ACCRETED STARS 

In order to examine the accretion history of galaxies, we re- 
strict our sample to galaxies with robust merger histories. 
We require that galaxies have a stellar mass M* > 1O 1O M , 
and at least 100 dark matter particles. We also require either 
that a galaxy can be traced back to z > 2.0 or that its small- 



© 2011 RAS, MNRAS 00Q.ITJfT8l 



Accretion vs. in- situ star formation 7 



est progenitor of the galaxy has a stellar mass M* < 1O 9 M . 
This removes galaxies which have gaps in their merger his- 
tory, but retains galaxies which have grown significantly 
since z — 2. The sample includes 443 galaxies from the clus- 
ter box and 168 galaxies from the void box. These repre- 
sent ~ 85% of the z — galaxies which satisfy the stellar 
and dark matter mass constraints. We test that this smaller 
sample is not a biased subset of the larger sample using the 
Kolmogorov-Smirnov (KS) test. We find the probabilities 
that the distributions of stellar mass, stellar age, and ac- 
creted fraction are the same for both samples are 1.00, 0.99, 
and 0.98, respectively. 

The top panel in Figure [4] shows the fraction of the final 
stellar mass which was accreted as a function of stellar mass. 
This figure (and all subsequent figures) only include galaxies 
for which we have robust merger histories. We find that the 
accreted fraction increases as a function of stellar mass from 
17 ± 7% for 10 10 < M* < 10 n M o to - 40 ± 15% for the 
most massive galaxies in the cluster box (M* > 1O 12 M ). 
This value is in agreement with estimates of the merged 
stellar mass for massi ve cluster galaxies from observations 
of paired galaxies (|Lin et alj|2010h . Although the sign of this 
trend is consistent with that found in OslO, both the slope 
and the zero-point are significantly smaller in our sample. 
The causes of this discrepancy will be discussed further in 
^3.11 We find no trend with the dynamical age of system and 
the ac crete d fraction, in c o ntrast to the results o f lFont et al. 
(|201lh and IZolotov et al.1 (|2009h . The excess of late in-situ 
star formation in our simulation may modify the dynamical 
ages of the galaxies and alter any trend in accreted fraction 
with age. 

The trend of accreted fraction with stellar mass is sen- 
sitive to the excess of massive galaxies in the simulations. 
If we restrict the sample to galaxies with M* < 10 11 , the 
slope of accreted fraction with mass is consistent with zero. 
Furthermore, overmerging will modify the slope of this re- 
lation. If the real relation between log M* and the accreted 
mass fraction, has a small or negative slope, then additional 
merging will tend to increase the slope; if the real slope is 
very steep, then overmerging will tend to flatten the slope. 
The value for which the slope remains unchanged is ~ 0.09, 
and depends only weakly on the stellar mass and average 
merger ratio. This value is close to the slope we measure 
slope, and suggests the slope in Figure [4] may be an artifact 
of overmerging, and the real slope could be either shallower 
or steeper. 

Direct comparisons of the merged fraction with obser- 
vations is difficult. Although there are man y studies of the 
merg e r rate as a function of redshift (e . g. iLe Fevre et al. 
2000; iBundv et all 1 20041; [Lin et al l liooi; iKartalteoe et al 
20071: lLotz et al.ll2QQ8l : Ijpgee et alJl2009l : iBridge et alJl20lA 
see ILotz et al.l 12011 and references therein), these stud- 
ies suffer from many observational biases, including diffi- 
culty in selecting comparable descendant /progenitor pop- 
ulations. Nonetheless, studies of paired galaxies indicate a 
weak positive tre nd in the major merger rate with pro gen- 
itor galaxy mass (IBundv et all 120091 : IParg et al ] l20ld . but 
see IBridge et alJl20ld ). This is compatible with our results 
of an increase in the accreted fraction with stellar mass, 
provided that the average merger ratio does not decrease 
precipitously with galaxy mass. 



3.1 Comparison to Oser, et al. 2010 

The upper and lower panels in Figure [H show the accreted 
fraction of stars for the galaxies in this work and in OslO, 
respectively. Both studies find the same sign for the rela- 
tionship between accreted fraction and total stellar mass, 
but OslO find an accreted fraction about twice as large as 
that found in this work. Furthermore, the slope of the OslO 
relation is about twice as steep (although the mass range of 
their sample is smaller). We suggest that these differences 
in accreted fraction can be explained by several major dif- 
ferences in the simulations in the resolutions, the feedback 
models, and the star formation efficiencies. 

As shown in Figure [5l even the most massive galaxies 
in our work experience significant late star formation. The 
median z — SFR for galaxies with stellar masses above 
10 11 Mq is 5.2 MQyr -1 , above the observed star formation 
rate for most present-day massive galaxies. This ongoing 
star formatio n contributes sign ificantly to the in-situ frac- 
tion of stars. iNaab et al.1 (2007) have shown that higher res- 
olution can decrease the amount of late time star formation 
(see their Figure 2). For the lower resolution simulations, 
the z — SFR is ^ 5 Moyr -1 . For smoothed particle hy- 
drodynamics (SPH ) codes, like the GADGET code used in 
Naab et a 1. (2007) and OslO, the l ow resolutio n suppresses 
the Kelvin-Helmholtz instability ([Agertz et alj|2007n , which 
breaks up clumps of cold gas before they sink to the cen- 
ter of a galaxy and form in-situ stars. At higher resolutions 
(ra* < 10 5 Mq), the final star formation rates is less than 
1 Mp/yr, even without AGN or Type I supernovae feedback 
([Johansson et al.ll2012h . 

Although the AMR code used in this work does not suf- 
fer from the same problem as SPH codes at low resolution, 
low resolution AMR simulations will suppress fragmenta- 
tion of infalling gas clouds. These gas clouds will survive 
to the center of the nascent galaxy, where they will form 
stars in-situ at late times. The high resolution simulations 
of OslO use particle masses of gas = 6.0 x 1O 6 M and 
^dark = 3.6 x 1 7 Mq. Although n ot the highest mass res- 
olution tested in lNaab et al.l (|2007h , this particle resolution 
is large enough to suppress most spurious late star forma- 
tion. Furthermore, the simulation in OslO have a factor 3 
more mass resolution in dark matter than the simulations 
examined in this work. Thus, we expect additional late in- 
situ star formation in our simulations, and, therefore, a lower 
accr e ted fr action when compared to the simulations of OslO. 
ICe3 (|2012T ) have performed spatial resolution tests for the 
cluster box and find that the overall convergence is quite 
good. However, they find that the metallicity in damped 
Lyman- alpha systems increases by around 0.2 dex for an 
AMR simulation with spatial resolution better by a factor 
of two. Therefore, we expect that improved spatial resolu- 
tion will alter the star formation and feedback processes, 
changing the the in-situ and accreted fractions. The same is 
certainly true for improved mass resolution, which will be 
tested in future work. 

Both the simulations presented here and those in OslO 
suffer from an overabundance of star formation, and a deficit 
of feedback. The differences in the feedback and cooling 
models may help explain th e difference in accreted fr action 
between the two simulations. IScannapieco et al.l (|201lh show 
that the addition of metal-line cooling to GADGET3 (used by 
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OslO) can help decrease the final stellar mass of galaxies by 
~ 30%, but does little to the star formation history. There- 
fore, we expect additional cooling will h ave little effect on 
the fr action of accreted stars. However, Scannapieco et all 
show that different feedback models can greatly af- 
fect the star formation history in simulations. For vari- 
ous models of fee d back and cooling in SPH simulations, 
IScannapieco et aD (|201lh find that the redshift by which 
half of the stars have formed ranges from ^ 1.5 — 4.0, com- 
parable to the redshift range of the median stellar age from 
OslO. Additional supernovae feedback, especially the addi- 
tion of supernovae winds, tends to decrease the median red- 
shift of star formation. We expect that decreased early star 
formation will decrease the accreted fraction, as fewer small 
systems will form stars at high redshift to be accreted later. 

In the work presented here, star formation occurs much 
later than in the work of OslO; Figure [5] shows half the 
stars in our simulations are only formed by z ~ 1, in agree- 
ment with results from other AMR codes (Sca nnapieco et aD 
l201lh . As mentioned above, late star formation will lead 
to an increased in-situ fraction in massive galaxies. We ex- 
pect that additional feedback will decrease star formation 
at late times and shift the median redshift of star formation 
to higher z. This will tend to decrease the in-situ fraction. 
Thus, improved feedback should bring our results and those 
of OslO into better agreement. 

In addition to differences in resolution and feedback, dif- 
ferences in star formation efficiency between the two sets of 
simulations affect the fraction in-situ versus accreted stars. 
iNaab et"al] ([2007) have shown that the total stellar mass is 
not sensitive to the star formation timescale. However, in 
simulations with short star formation timescales, stars will 
form in small clumps outside the galaxy and will be later 
accreted, while in galaxies with longer star formation times, 
small clumps will instead be accreted as gas and later form 
stars in-situ. Thus, the total number of baryons accreted by 
the galaxy remains the same, but whether the baryons are 
accreted as stars or gas depends strongly on the star for- 
mation timescale. This effect is at least comparable to the 
effects from feedback on the accreted fraction. In comparing 
two grid-based codes with different star formation efficien- 
cies, the redshift by which half the stars have formed moves 
from z ~ 2 for an efficiency o f 5% to z ~ 1 for an efficiency of 
1% (|Scannapieco et al.ll201lh . comparable to the change in 
the median redshift of star formation from different feedback 
models. 

Ignoring the return of mass to the gas phase from high- 
mass stars, the star formation rate can be written as 



dt 



n M gas 
£dyn 



(2) 



where tdyn is the dynamical time of the gas forming stars 
(see equation [T]) . Thus, 1/C* x tdyn is the star formation 
timescale, tsF- Simulations with large values for C* have 
short tsF and therefore make stars efficiently. In the simula- 
tions examined in this work, C* — 0.03, in agreem ent with 
theoretical calculations (Krumholz & McKee 2005) and ob- 
servations (|Kennicutt| ["l998; Evans et al. 2009). In OslO, the 
value for C* is 0.083, which means for a given density gas 
clump, star formation is ~ 2.8 times faster in the simula- 
tions used by OslO. This difference in star formation rate 



becomes especially important for small clumps of baryons 
being accreted onto larger galaxies. 

We can compare the star formation time scale to the 
timescale for material to fall in from the virial radius. The 
infall time is given by the dynamical time of the dark matter 
halo, tviriai oc 1/ V^Pviriai = 1 / ^OOGpcrit • Thus, the ratio 
of the star formation time to the infall time is given by 



tsFy/tvirial 



V200Gp C rit 

gas 

1 /300 



(3) 



-^0 tdy 



0.05/i 



tdyn 



100 Myr 



For tdyn = 100 Myr, the ratio of the star formation timescale 
to the infall timescale is 0.6 h and 1.7 h for C* = 0.083 and 
0.03, respectively. In the first case, an infalling clump of 
gas above the star formation threshold density will convert 
~ 91% of its gas to stars, while in the second case only 
57% of the gas will be converted into stars, and the accreted 
fraction will decrease by about 40%. The difference in C* 
between the two simulations partially explains the differ- 
ence of a factor of 2 between the accreted fractions found in 
this work and OslO. It also highlights the importance of the 
star formation efficiency in studying the merger histories of 
galaxies. Unlike dark matter halos, which always undergo 
'dry' mergers, baryons can merge either as stars or gas. The 
relative fraction of stars to gas can have profound effects on 
the observable properties of galaxies. 



3.2 Properties of accreted and in-situ stars 

In addition to examining the total contribution of in-situ 
and accreted stars, we can also examine the mean proper- 
ties of the accreted and in-situ fractions as a function of final 
galaxy stellar mass. The top panels in Figures [2] and [3] show 
the mass-weighted stellar ages and stellar metallicities for 
the in-situ and accreted stars as a function of z — galaxy 
stellar mass. It is immediately evident that the accreted stars 
are older (by 2 — 2.75 Gyrs) and more metal-poor (by ^0.15 
dex) than the stars formed in-situ. Figure [5] shows distribu- 
tions of stellar ages for galaxies divided into three stellar 
mass bins. The vertical lines show the median ages for the 
in-situ, accreted, and total stellar masses. Again, it is clearly 
evident from these figures that the accreted stars are older 
than the in-situ stars. The plots in this figure as well as Fig- 
ures [6]— [9] weight all galaxies in the cluster and void boxes 
equally. If we weight th e galaxies a ccording to the fit to the 
star formation history (Ceni l2011bh , only the results for the 
lowest mass bin change, and the differences in the median 
values are very small. Furthermore, by dividing the galaxies 
into mass bins, we see in Figure [5] that there is more late- 
time star formation in lower mass (M* < 10 11 M ) than 
i n high er mass galaxies, consistent with the results of ICenl 
()2011br ) using the same simulation presented here. Exam- 
ining the trends in median stellar age for the accreted and 
in-situ fractions as a function of mass reproduces the trends 
shown in Figure [2] namely, while the in-situ and accreted 
stellar ages only change slightly with stellar mass, the me- 
dian age of all the stars increases since the accreted fraction 
grows with stellar mass. Finally, the difference in stellar age 
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Figure 5. The probability distribution of stellar ages for galaxies 
grouped by mass. The black solid line is the total stellar age 
distribution. The blue dashed line is the accreted stars and the 
red dash-dotted line is the in-situ stars. The similar vertical lines 
show the median age of the stars in each category. For all mass 
bins, the accreted stars are ~ 2 Gyrs older than the in-situ stars. 

between the accreted and in-situ stars is roughly indepen- 
dent of total stellar mass. If we exclude galaxies more mas- 
sive than 10 111 M© and only examine the top panel, the 
accreted stars are still more than 2 Gyrs older on average 
than stars formed in-situ. 

Since the age distributions of the in-situ and accreted 
stars are different, it is not surprising that the metallic- 
ity distributions are different as well (see Figure [6]). For 
the most massive galaxies, the accreted stars are 0.12 dex 
more metal-poor than the stars formed in-situ. The sign of 



Figure 6. The probability distribution of stellar metallicities for 
galaxies grouped by mass. The lines are as in Figure [5] At all 
masses, the accreted stars are more metal poor than in-situ stars. 

this difference is in agreement with simulation results from 
iFont et al.l (|201ll ) for disk galaxy halos, but the metallicity 
difference found in this work is smaller. However, as noted 
above, the stellar metallicities in the simulation are sensitive 
to the supernovae feedback prescription, and are highly un- 
certain. Nonetheless, it is unlikely that improved feedback 
would reverse the sign of the difference in metallicity be- 
tween accreted and in-situ stars. Additionally, the difference 
in metallicity between accreted and in-situ stars is largest 
for the lowest stellar mass bin, where the mass function of 
galaxies agrees with observations. However, additional su- 
pernovae feedback may play an important role here, as well. 
The bimodal distribution of the stellar metallicities in 
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the low mass galaxies (see the top panel of Figure [6]) may 
be partially due to over- weighting the cluster galaxies rela- 
tive t o the void galaxies. If we weight the galaxies as in Cen 
(|2011bh (the void galaxy contribution is increased by a fac- 
tor of 6) , the metallicity distribution in the top panel of Fig- 
ure [6] is still bimodal, but the low metallicity peak is nearly 
equal in height to the high metallicity peak. This remaining 
bimodality may be due to enhanced late star formation in 
metal-rich gas, which is not observed in present-day galax- 
ies, and probably due to inadequate resolution and stellar 
feedback. 

In addition to differences in the median stellar metallici- 
ties for accreted and in-situ stars, the dispersions are also dif- 
ferent; the accreted stars have a larger metallicity dispersion. 
For the three panels in Figure \6\ the inter-quartile ranges of 
the in-situ stellar Z/Zq are 0.26, 0.28, and 0.27 dex. The 
same ranges for the accreted stars are 0.60, 0.41, and 0.34 
dex. The differences in metallicity distribution width are due 
to the inhomogeneous make-up of the accreted material. The 
accreted stars are formed in a variety of systems with differ- 
ent metallicities, while the in-situ stars are formed in a single 
stellar system. The same difference in me tallicity dispersion 
is seen in the Galactic disk and halo fe.g. lBeers et al.| [l985) 
and is used a s evidence for the accretion origin of the G alac- 
tic halo (e.g. Searle &; Zin nl ll978l : lUnavane et"Hlll996h . 

We also examine the spatial distribution of in-situ and 
accreted stars in galaxies at z = 0. Figure [71 shows the 3- 
dimensional average stellar mass profiles for the accreted 
and in-situ stars in three differ ent mass bins. A s expected 
(and as shown in the work of lOser et al ] <|2012h ), the half- 
mass radius of the accreted stars is larger than the half-mass 
radius of the in-situ stars by approximately a factor of 2 for 
all three mass bins. Some of the accreted material on the 
outskirts of galaxies may be due to our choice of density 
threshold in the HOP grouping algorithm; a lower threshold 
will lead to larger galaxies with larger accreted fractions. 
In order the check how much this affects our results, we 
show the median radius of the stars accreted from resolved 
galaxies only in Figure This radius is still larger than the 
median radius of the in-situ stars. Furthermore, this radius 
represents the minimum half-mass radius for the accreted 
material, as additional accretion from unresolved galaxies 
will tend to increase this radius. 

Trends with accreted fraction and radius suggest that 
the average merger event is not an equal mass- merger, 
after which the stars would be well-mixed, but rather a 
minor merger in which low density material is added to 
the outside of a massive galaxy. It also suggests that mi- 
nor mergers can substantially increase the half-light radius 
of massive galaxies, leading to the growth of early- type 
galaxies at late times. This is in agreement with both ob- 
serva t ions teezanson et alj [20091; Ivan Dokkum &; Brammerl 
l20ld : iNewman et alJl20ia iTal et alJl2012h and other sim- 
ulations (iGallagher fe Ostrikerll 19721: iBovlan-Kolchin fe Mai 



l2004l ; iNaab et al.ll2007l . 120091 : lOser et aljl2012h . It also fol- 
lows from simp le, semi-quantitative arguments based o n the 
virial theorem (|Naab et al.ll200d : iBezanson et alJl20Q9h . Be- 
cause the accreted material has lower metallicity than the 
in-situ material, we expect galaxies to have negative metal- 
licity and colour gradients. In projection, the average metal- 
licity gradient for the galaxies with robust merger histories is 
dlog(Z/Z )/dlog(i?/i? e ff) = -0.38 ±0.23, where i? e ff is the 



projected half-mass radius. Furthermore, the differences in 
the distributions of stellar metallicity between the accreted 
and in-situ stars lead to an increase in the metallicity dis- 
persion as a function of radius. The standard deviation of 
log Z/Zq has an average slope of da\ og z/z Q /d \og(R/R e ft) — 
0.28 =b 0.18. There are many mechanisms besides accretion 
that will contribute to an increase in metallicity dispersion. 
Indeed, we find the in-situ stars alone also have an increase 
in metallicity dispersion with radius only slightly smaller 
than the increase reported above. 

Observations of these gradients are difficult be- 
cause the low surface brightness in the outskirts of 
galaxies. There are numerous me asurements of colour 
gradients in early-type galaxies (Ide Vaucouleursl 
lEisenhardt et al.ll2007l: iTortora et alJl20ld: ISuh et al. 



1961; 



201C 



|Tal fe van Dokkumll2oTTI;lLa Barbera et al.ll201ll : lGuo et all 
1 20 111 ; Gonzalez- Perez et alj l201lh , although only a few 
studies extend beyond the half-light radius of the galaxy. 
Furthermore, colour measurements cannot be used to 
distinguish between metallicity and age gradients. More 
recently, studies of the metallicity gradients of globu- 
lar cluster systems in early-type galax ies support theo- 
ries of two-phase galaxy formation ([Arnold et al.l 1 20 111 ; 
iForbes et alJuOlU ). Direct measurements of metallicity gra- 
dients in the stellar populations of early type galaxies have 
been done using long slit spectroscopy (Davies et al. 1993; 



„ ___ _ £ ,-jpy ^ _ 

iMehlert et al.ll2003l : lBrough et al.ll2007l:ISpolaor et alJ201Ch. 
and integral-field-un its jRawle et all 2008: Kuntschner et al. 
2010; Greene et al. l l20ia T Only a few measurements focus 
on gradients beyond the half-light radius (iFaber et al] 19771 ; 
IFoster et all l2Q09l : ISoolaor et atll20ld : iGreene et alJl2012h . 
These studies find metallicity gradients (but not age gradi- 
ents) which indicate that accretion of low metall icity systems 



is important in the build-up of massive galaxie s (IFoster et al 
2009; So olaor et al .11 20 id : IGreene et alJl2oT3 ). IGreene et al 



(|2012h find an average slope of [Fe/H] « -0.1 from 0-2.5 
i? e ff for 8 massive early type galaxies, in good agreement 
with the average slope we measure in the simulations. Us- 
ing the iron and magnesium abundances in the outskirts of 
galaxies, they obtain an average mass ratio for the merger 
events of 10 : 1, which agrees roughly with our estimate for 
the mean merger ratio below (see . 

Although the accreted material is made up of older 
stars, it is not added to the galaxies until late times. Fig- 
ure [8] shows the mass- weighted average mass assembly for 
galaxies in three mass bins. For the central mass bin, half 
of the accreted stars are formed before z ~ 1.6, but it takes 
until z ~ 0.6 for half of these stars to actually be accreted. 
Furthermore, Figure [8] shows that the in-situ stars are in 
place before the accreted material is added. The difference 
in assembly times between the accreted and in-situ stars is 
largest for the lower mass galaxies. This i s consistent with 
the picture of two phase galaxy formation ([Oser et al.ll20ld ; 
INaab et ai]l2007h : galaxies first undergo a phase of intense, 
in-situ star formation, followed by a phase of accretion of 
old, less massive, and therefore more metal-poor systems. 

3.3 Mergers and environment 

We can also use the galaxies from the cosmological simu- 
lations to examine the properties of mergers as a function 
of environment. This analysis is not possible for simulations 
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Figure 7. The cumulative, 3-dimensionally averaged radial stel- 
lar mass profiles for the galaxies in three mass bins. The lines are 
the same as in Figure [5] The vertical lines show the 3-dimensional 
half-mass radii for the in-situ, accreted and total stellar mass 
profiles. The green dotted lines show the half-mass radius for ma- 
terial accreted from resolved galaxies only. These represent the 
minimum radius for the accreted stars. The accreted stars are 
consistently found at larger radii than the stars formed in-situ. 



which only include individual halos. Studies of paired galax- 
ies have f ound a dependence on the merger rate with envi- 
ronment rtMcIntosh et all 120081 : iDarg et all 12010k iLin et al.1 
120101 : iTonnesen fe Cenll201lh . This is expected, as higher 
density environments have more galaxies and higher poten- 
tial for merging. However, part of this trend may be due 
to the correlation of environment and galaxy mass, which 
is related to the accreted fraction. By selecting a sample 



Figure 8. The stellar mass assembly histories for galaxies in three 
mass bins. The lines are the same as in Fig. \E\ The vertical lines 
show the time at which 50% of the in-situ, accreted, and total 
mass is in place. While the accreted stars are typically formed at 
z ~ 1.5, they are assembled at z ~ 0.6. 

of galaxies of similar mass, we can eliminate trends in the 
accreted fraction with mass and isolate the trends with en- 
vironment. 

We define environment as the 3-dimensional 
galaxy density to the 5 th nearest neighbor (p 5 = 
6/(distance to 5 th neighbor) 3 [/i 3 Mpc -3 ]). The densi- 
ties used are an average of the z — and z — 0.05 measured 
densities in order to eliminate some of the noise introduced 
by the nearest neighbor density measure. We have restricted 
the sample to 109 galaxies with stellar masses between 
6 x 10 10 M < M* < 3 x 10 11 M . Overall, we find that 
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the accreted fraction is not a strong function of present-day 
environment; the average accreted fraction increases from 
~ 0.20 for field galaxies to ~ 0.2 5 for cluster galaxies . 
This contradicts the observations of iMcIntosh et all (|2QQ8h , 
which show an increase in the merger rate for group 
galaxies by a factor of 2-9 over the rate for field galaxies. 
The difference is due in part to our decision to study trends 
with environment at const ant stellar mass , elim inating the 
larger trend. Furthermore, [Mcintosh et al.l (2008) use group 
and cluster membership as a measure of environment. This 
measure is more sensitive to the near-field environment 
of a galaxy than ps, and therefore, it is probably more 
closely correlated with the number of mergers a galaxy has 
undergone. 

Figure [9] shows the stellar age distributions for galaxies 
in three density bins, which correspond to the local density 
for field, group, and cluster galaxies. Typical separations 
between galaxies in each environment bin are 5 /i -1 Mpc, 
1 /i -1 Mpc, and 0.5 /i -1 Mpc, respectively. These bins were 
chosen to have equal numbers of galaxies in the mass range 
6 x 10 10 M o < M* < 3 x 10 n M o . However, the results are 
not sensitive to the exact locations of the bin separations. 
As expected, the star formation histories of the galaxies in 
low density environments are far more extended. For field 
galaxies, both the in-situ and accreted stars are younger 
than in group and cluster galaxies by ~ 2 Gyr. In fact, not 
only does the in-situ star formation occur later in low density 
regions, but the accretion also occurs ~ 1.5 Gyr later than 
in the higher density regions. In every respect, low density 
regions evolve more slowly than high density regions. 

Figure [TO] shows the total star formation histories for 
galaxies in low, medium, and high density regions. These 
histories are normalized to the local dark matter mass in 
each region, which is accomplished by weighting each galaxy 
by the dark matter mass enclosed in a sphere centered on the 
galaxy and extending to the galaxy's fifth nearest neighbor. 
As expected, star formation in galaxies in the lowest den- 
sity environments dominates today, while, at high redshift, 
star formation occurred mainly in the mid- and high density 
regions. Thi s is c o nsisten t with the star formation histories 
presented in lCenl (|2011bh . 



4 ACCRETED STARS AND MERGERS 

The accreted fraction of the final stellar mass is added to 
the galaxies through mergers. Using the simulations we can 
determine the rate at which mergers occur and the distri- 
bution of mass ratios of these mergers. Since the accreted 
fraction in most of our galaxies is less than 30% we expect 
that minor mergers (merger mass ratios less than 4:1) will 
dominate the accretion history of galaxies in our sample. 
In the following sections, we compute the merger rate for 
our sample as a function of redshift and the average merger 
ratios for the most massive galaxies in the simulation. 

4.1 Merger rates 

The slopes of the in-situ and accreted mass assembly histo- 
ries in Figure[8]are the star formation rates and merger rates, 
respectively. The star formation rate for the massive galax- 
ies (logM*/M > 12.2) peaks at z ~ 2.5, when the star 



-3.0<logp 5 <-0.6, N = 37 
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Figure 9. The probability distribution of stellar ages for galaxies 
grouped by environment (ps)- The galaxies all have stellar masses 
between 6 x 1O 1O M < M* < 3 x 1O 1O M . The lines are the 
same as in Figure [5] Although the accreted fraction is only a 
weak function of environment, the median stellar age of the in- 
situ stars is a strong function of environment. 

formation rate is high enough to double the galaxy mass 
within a Gyr. The accretion rate for the massive galaxies 
peaks later, at z « 1, with an average growth rate from 
accretion of ~ 15% Gyr -1 . At z = 0.25, the galaxies with 
M* > 3 x lO n M are adding stellar mass in mergers at a 
rate of ~ 2.6% Gyr -1 . This includes growth from unresolved 
mergers in the simulations. 

The simulation merger rate is approximately a fac- 
tor of two larger than the observed merger rate for lumi- 
nous red galaxies (LRGs) of 1.7%h Gyr -1 , where Hq — 
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Figure 10. The star formation rates in different environments, 
for galaxies with stellar masses 6 x 1O 1O M < M* < 3 x lO n M 
at z = 0. The different curves are normalized to the average dark 
matter mass associated with each galaxy, as given by 1/5 of the 
dark matter mass to the fifth nearest neighbor. At early times, 
the star formation occurs in the densest regions, while at late 
times the field galaxies dominate the star formation. 



.20091 : 



and quantified morphological disturbances (e.g.! ^otz et aD 
120081 : IConselice et all [20091 : lLopez-Saniuan et alJl2009h , re- 
spectively. 

We can also compare the evolution of the merger rate 
as a function of redshift to the observations of merger rate 
evolution. However, care must be taken to ensure the selec- 
tion of galaxies from the simulations and observations are 
similar. lLotz et aD ([201 ll ) argue that galaxies selected based 
on luminosity or mass cuts do not give the same evolution as 
samples selected based on constant galaxy number density, 
as the progenitor and descendant galaxies are different in the 
mass and luminosity selected samples. In our simulations, 
the merger rate per galaxy per Gyr for galaxies with stellar 
masses above 10 10 M grows as (1 + z ) 2 - 0±0 - 25 for z < 1.5. 
Thus, we expect the merger rate to be an incre asing function 
of re dshift, i n agreement with observations ([Bridge et al.l 
l20ld . but see iBundv et alJ l2009h . The slope of the galaxy 
merger rate with redshift is similar t o the slope of the dark 
matter merger rate w i th redshift (e.g. iGottlober et all l200ll : 
iFakhouri fe Mai 120081 : iFakhouri et~afl l2010h . However, the 
simulations clearly show systems in which the dark matter 
halos merge well before the galaxies merge. Cosmological 
hydrodynamical simulations, such as those used here, can 
be used to measure the time lag between dark matter halo 
mergers and galaxies mergers, which we plan to pursue in 
future work. 



100/i km s" 1 Mpc" 1 , found bv Masiedi et al.l <j20Q8h using 
the small scale cross correlation function to determine the 
close pair fraction in SDSS. This observed rate includes 
growth from minor mergers since the maximum luminos- 
ity ratio between th e LRGs and other g alaxies is ~ 50:1 in 
SDSS observations ([Masiedi et al.|[200a i. Since the massive 
galaxies in the simulation suffer from overmerging, we ex- 
pect the simulation merger rate to be inflated compared to 
the real value. 

Comparing our simulated merger rates to other 
observed merger rates based on paired galaxies or 
morphologically-disturbed galaxies is not straightforward; 
it requires knowing the average merger mass ratio and the 
merger rate of a single population of galaxies as a function of 
time. Assuming an average merger mass ratio between 0.15 
and 0.3 (compatible with our results from fc|4. 3[) , the merger 
rate as meas ured by qu antified morphological disturbances 
(G-M 2 o) in lLotz et al.l (|201lh is 7.5-15% Gyr -1 at z = 0.2 
for galaxies at number densities above - 6 x 10" 3 Mpc -3 . 
This is considerably larger than both the me rger rate mea- 
sured in our simulations and that measured bv lMasiedi et al.l 
(2008), suggesting that the observational merger rate is still 
uncertain. Using h ighe r resolu tion simulations with the same 
initial conditions JCenl ([201 ldh find that the observationally- 
deflned merger rate from close pairs of galaxies does ex- 
ceed the physical merger rate for simulated galaxies. This 
over-estimate is at least partially correcte d for by includin g 
an observable timescale for the merger ([Lotz et al.ll201lh . 
Despite this correction, the observed merger rates are still 
large. The total merger rate from our simulations is com- 
parable to the observe d major merger rates from studies of 
close galaxy pairs (e.g. IBundv et al.l 2004; Ka rtaltepe et al.l 
120071 : Ide Ravel et al. 2009; Bu ndv et al] 120091 ): visible dis- 



4.2 Expected merger ratio 

In addition to studying the rate of merging, the simulations 
also allow us to examine the average merger mass ratios as 
a function of galaxy mass. We define the merger ratio, /x, as 
the mass of the accreted galaxy divided by the mass of the 
progenitor galaxy. By definition, fi is always less than one. 
We first calculate the expected merger ratio given the galaxy 
mass function, and compare this expected value with the 
value from simulations. Similar calculations for the limiting 
case o f small merger ratios are done in OslO and lNaab et al.l 

(Hooi). 

Ignoring effects of dynamical friction, gravitational fo- 
cusing, and local density, the mean merger mass ratio for a 
galaxy is the average mass of less massive galaxies divided 
by the parent galaxy mass. This can be calculated from the 
mass function. Given a mass function $(M)dM, the prob- 
ability of a galaxy with mass Mo accreting a galaxy with 
mass in the range (M a , M a + dM) is 



P(M a)d M = J^ dM 



(4) 



and the number-weighted and mass-weighted average 
merger mass ratios are simply 

(fi n ) — I mP(m)dm/ I P(m)dm , and 

J M min J M min 

(fim) — / m 2 P(m)dm/ I mP(m)dm , (5) 

J M min J M min 



respectively. For a Schechter function, P(M) oc 
(M/M*) a exp~ M/M *, and the number-weighted and 
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Figure 11. The mean merger ratios as a function of stellar 
mass. The red points are from the cluster box and the blue 
crosses are from the void box. The top panel shows the number- 
weighted mean merger ratio, while the bottom panel shows the 
mass-weighted mean merger ratio. The solid line is the expected 
mean merger ratio computed using the unweighted mass function 
from the simulation. The dotted line shows the expected mean 
merger ratio for a Schechter function galaxy mass distribution 
witha = -1.16 and log M* = 10.8. 



mass-weighted average merger ratios can be evaluated as 

(/in) = 1/Mo x (r (a + 2, M min ) - T (a + 2, M )) / 

(r (a + 1, M min ) - T (a + 1, Mo)) , 
(Hm) = 1/Mo x (r (a + 3, M min ) - T (a + 3, M )) / 

(r(a + 2,M min ) -r(a + 2,M )) , (6) 

where masses are given in units of M* , and we have included 
a lower limit on the mass function, M m i n . The dashed lines 
in Figure ITT1 show the numerical values for (/i n )and (/i m )for 
a = —1.16 and logM*/M = 10.8, the Schechter fun ction 
parameters fit to the SDSS sample (|Li fc White! 1 20091 ). The 
integrals are stopped at and M m in = 10 9 M . This pro- 
cedure yields simil ar results to the exp ected merger mass 
ratio calculated by I Hopkins et al.l (|2010h using N-body sim- 
ulations for the merger trees and se mi-empirical methods 
for determining the galaxy properties. iHopkins et al.l (|2010l ) 
find that (/i m ) peaks around M* w 10 11 M , and decreases 
for both larger and smaller galaxy masses. 



4.3 Mean merger mass ratio 

In fc|2.2l we show that the simulated mass function devi- 
ates from the observed mass function below M* < 1O 1O M0. 
Therefore, in this section, we only examine mergers between 
galaxies with stellar masses greater than 10 10 M© at z < 2.5. 
We require that the final galaxy mass is larger than 10 11 Mq , 
thus allowing us to resolve mergers with a mass ratio of 
fi < 0.1 at z = 0. This limits the sample to the most mas- 
sive galaxies which suffer from over-cooling and too much 
star formation, making comparisons to observations inef- 
fective. This sample includes 85 galaxies from the cluster 
box and 15 galaxies from the void box. Figure [11] shows 
the mean mass-weighted and number-weighted merger ra- 
tios as a function of z — stellar mass. The large triangles 
denote median merger ratios and do not show a significant 
trend with galaxy mass, although the maximum merger ra- 
tio increases with decreasing stellar mass. The averages are 
consistent with the expected merger ratio from the mass 
function in the simulation (solid line). The excess of mas- 
sive galaxies due to over-cooling makes the expected merger 
ratios in the simulation deviate from the merger ratios pre- 
dicted by a Schechter function galaxy mass distribution (cf. 
Figure [1]). The 3-a clipped average values for (/i n ) and (/i m ) 
are 0.20 =b 0.16 and 0.26 =b 0.21, respectively. These results 
are roughly consisten t with the merger ratios calculated in 
Hopkins e t all (|20ld) using semi-empirical techniques. 

lOser et all (2012) examine the average merger ratio for 
the same set of galaxies used in OslO. They also find little 
trend with merger ratios and stellar mass, but their aver- 
age merger ratios are lower than the ones we find here and 
closer to the analytically derived values. Because the simu- 
lations in OslO are of higher mass resolution, they are able 
to resolve mergers down to lower mass ratios. Therefore, it 
is not surprising that they find a substantially lower value 
for (/i n ), which is affected by many low mass mergers, while 
their value for (/i m ) is comparable to the (/i m ) we mea- 
sure. Because we are restricted to mergers with mass ratios 
greater than 10:1, the plots in Figure [TT] only account for 
75% of the accretion onto the galaxies shown. The other 
25% of the accretion is from galaxies below our mass cutoff 
and other unresolved accretion. Accounting for this accre- 
tion sets a lower limit on (fa) of 0.20, in agreement with 
the results from lOser et al.l (l2012h. Fin ally, the low value 
for (/i n ) reported by Oser et al.l (|2012h is consistent with 



the differences in star formation efficiency in the two works; 
the h igher star formation efficiency in OslO and lOser et al l 
f|2Q12h should yield more small stellar systems than in our 
simulations, which will decrease (/in). 



5 DISCUSSION AND SUMMARY 

The two-phase picture of galaxy formation put forward by 
OslO provides a useful framework for studies of galaxy evo- 
lution. We apply this framework to study the galaxy merger 
histories of more than 600 simulated galaxies with M* > 
10 10 M©. In this work, we corroborate the two-phase model 
for galaxy formation. The first phase consists of a period of 
in-situ star formation, which occurs around 0.75 < z < 1. 
This in-situ phase accounts for 60-90% of the star forma- 
tion in galaxies with stellar masses above 10 10 M . The 
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Figure 12. Distributions of stellar ages and host galaxy masses for accreted (blue contours) and in-situ (red contours) stars. The x-axis 
shows the mass of the galaxy in which a given star particle formed. The left panel shows the distribution for a galaxy which underwent 
a major merger at z ~ 0.4. The center panel shows a galaxy which had no major mergers. The right panel shows the distribution of 
stars for all massive galaxies in the simulation. The contours have been smoothed and enclose 20, 50, and 80 per cent of the stars. The 
differences between the accreted and in-situ stars are largest for galaxies with no major mergers. 



remainder of the stars are added in mergers. The peak of 
the merger activity occurs at z ~ 0.5, after the majority of 
the in-situ star formation has taken place. Our results are in 
good agreement with those from OslO and other simulations 
distinguishing between accreted and in-situ star formation 
in stellar halos, mainly around disk galaxies ( e.g. iFont et al.l 
l201ll ; lAbadi et al.ll2006l : IZolotov et aHl2009l l2010h . In par- 
ticular, the differences in stellar age, metallicity, and spa- 
tial distribution shown in Figures EH?] are in a greement with 
the di fferences found bv lFont et al.l (|2Q1 lh and lZolotov et all 

Although the accretion occurs at late times, the ac- 
creted systems are old. The accreted stars are on average 
2 Gyr older than the stars formed in-situ. This age differ- 
ence correlates with a median metallicity difference of 0.15 
dex. Since the accreted material comes from a variety of 
sources (i.e., a range of galaxy host masses), the distribu- 
tion of stellar metallicities is larger for the accreted stars 
than for the in-situ stars. However, the simulated galaxies 
do not conform to the mass-metallicity relation, suggesting 
that the stellar feedback model does not provide sufficient 
feedback, and that higher resolution is needed. Improved 
feedback prescriptions should lower the stellar metallicities 
for the less massive galaxies and the accreted systems, thus 
increasing the differences in metallicity between in-situ and 
accreted stars. This effect may be countered by a decrease 
in metallicity in the centers of galaxies, if higher resolution 
and better feedback prescription help shut off the excess late 
star formation currently seen in the simulations. 

Unlike in simulations, observations cannot easily distin- 
guish between stars formed in-situ and those added by ac- 
cretion. Nonetheless, we find indirect signs of the two phases 
of galaxy formation which are observable today. Because the 
accreted stars formed in small systems at early times, their 
stellar ages and metallicities are noticeably different from 



the stars formed in-situ. Furthermore, because the accreted 
stars preferentially reside in the outskirts of galaxies, we find 
strong gradients in the stellar populations of massive galax- 
ies which have been observed beyond the half light radius of 
early-type galaxies ([Foster et al.l 120091 : ISpolaor et al.l l2010l : 
iGreene et al1l2012h . 

As explained above, our results are consistent with mod- 
els of hierarchical formation of galaxies; large systems are 
assembled at late times from older, smaller building blocks. 
Figure [T2l is an illustration of the hierarchical growth. Here, 
we plot the stellar age against the stellar mass of the system 
(galaxy) in which each star was born, for the accreted and 
in-situ stars, separately. For a galaxy with no major merg- 
ers (central panel), the distribution of accreted and in-situ 
stars are well separated. On average, accreted stars form 
early in small systems, while in-situ star formation occurs 
later, when the host galaxy is more massive. On the other 
hand, the distributions of accreted and in-situ stars are very 
similar for a galaxy which underwent a major merger (left 
panel). This is unsurprising since two galaxies with roughly 
the same mass should have similar star formation histories. 

The overall trends in Figure [12] suggest that galaxy stel- 
lar mass and stellar age are anti-correlated. However, the low 
mass galaxies shown in this figure formed early and were 
subsequently accreted, at which point they stopped form- 
ing stars. These galaxies are not representative of present- 
day low mass galaxies, but are rather the building blocks 
of more massive galaxies. Indeed, the top panel of Figure [2] 
shows stellar mass and stellar age are correlated for galax- 
ies which survive to z = 0, in qualitative agreement with 
observations. 

Observations cannot directly reproduce Figure IT2l 
However, there may be a close mapping from the rela- 
tion between stellar ages and progenitor masses to the 
relation between metallicity and a-element enhancement 
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(jThomas et al.ll2005l : Ijohnston et al.ll2008T ). Old stellar pop- 
ulations have higher alpha-to-iron abundance ratios than 
younger stellar systems, while massive galaxies typically 
have higher metallicities than less massive systems. Figure 
IT2l shows that accreted stars come from older, less mas- 
sive systems, so we expect that accreted stars should be 
alpha-element enhanced compared to stars formed in-situ. 
This will yield radial gradients in abundance ratios in ad- 
dition to the gradients in metallicity discussed above. This 
stellar abundance space is used to map the accretion his- 
tory of the Galactic ha l o (IBl and-Hawthorn & Freeman 2003; 
Robertson et al. 2005; Font et al.l l2006b llal; Ijohnston et alj 



2008: IZolotoT et al. 201 j ). Usin gsimulations of a Milky Way 



analogue, iTissera et al . (2012) find that the accreted stars 
in the disc (15% of the disc stars) are alpha-enhanced and 
older than the stars formed in-situ, in broad agreement with 
our results shown in Figure [12] This type of analysis is being 
extended to galaxies besides the Milky Wa y by measuring a- 
to-Fe ratios as a function of galaxy radius (iKuntschner et all 
120101 : ISpolaor et all 120101 : IGreene et all l2012h . As detailed 
chemical and spatial information becomes available for the 
outer regions of more galaxies, it may be possible to differ- 
entiate between a minor merger and major merger driven 
galaxy formation history for individual galaxies. 



The results of this work agree broadly with the results 
of OslO based on a set of 40 galaxies simulated using GAD- 
GET, an SPH code, at higher mass resolution. We claim 
that most of the differences can be accounted for by the 
differences in star formation efficiency, feedback, and reso- 
lution between the simulations. The higher star formation 
efficiency in OslO yields a higher accreted fraction because 
small halos have time to form stars during infall from the 
virial radius of larger halos. These small stellar systems 
tend to decrease the number-weighted mean merger ratio 
(|Oser et al . 2012). In our simulations, the lower resolution 
tends to increase the in-situ star formation rate by prevent- 
ing gas clumps from fragmenting and forming stars before 
being accreted. These gas clumps survive to the centers of 
the galaxies, where they then form stars in-situ. The effects 
of resolution on the in-situ and accreted fractions will be 
checked in future AMR simulations run at higher resolu- 
tion. Together with the differences in feedback, these effects 
account for the difference of a factor of two between the 
accreted fraction observed in the SPH and AMR simula- 
tions. Despite these differences, there is an excellent quali- 
tative agreement between this work and the work in OslO. 
This shows that the two phase model for the formation of 
massive galaxies is robust and independent of the numeri- 
cal methods and subgrid physics. Therefore, the simulations 
can be used to predict observable consequences of the two 
phases of galaxy formation. In this work, we show that the 
differences in age, metallicity, and metallicity dispersion be- 
tween the accreted and in-situ stars yield radial gradients in 
these quantities for present-day massive galaxies. Currently, 
there is evid ence for radial trends in metallicity in ea rly-type 
galaxies fcf. IGreene et al.ll2012l ; ISpolaor et al.ll2010h , consis- 
tent with our predictions for two phase galaxy formation and 
further observational tests of this picture will be extremely 
useful. 
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